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ABSTRACT 

Context. Abundances of galaxies at redshifts z > 4 are difficult to obtain from damped Ly a (DLA) systems in the sightlines of quasars 
(QSOs) due to the Ly a forest blanketing and the low number of high-redshift quasars detected. Gamma-ray bursts (GRBs) with their 
higher luminosity are well suited to study galaxies out to the formation of the first stars at z > 10. Its large wavelength coverage 
makes the X-shooter spectrograph an excellent tool to study the interstellar medium (ISM) of high redshift galaxies, in particular if 
the redshift is not known beforehand. 

Aims. We want to determine the properties of a GRB host at z = 4.66723 from absorption lines. This is one of the highest redshifts 
where a detailed analysis with medium-resolution data is possible. Furthermore, we determine the dust extinction using the spectral 
energy distribution from X-rays to near infrared. Finally, we put the results in context to other GRBs at z > 4 and properties of (high 
redshift) QSO absorbers. 

Methods. The velocity components of the resonant and fine-structure absorption lines are fitted with Voigt-profiles and the metallicity 
determined from S, Si, Fe and O. Si IP together with the energy released in the UV restframe as determined from GROND photometric 
data gives us the distance of the absorbing material from the GRB. The extinction is determined from the spectral slope using X-ray 
spectral information and the flux calibrated X-shooter spectrum, cross calibrated with photometric data from GROND. We also collect 
information on all GRB hosts with z > 4. 

Results. We measure a relatively high metallicity of [M/H] = -1.0 + 0.1 from S, the distance of the material showing fine-structure 
lines is 0.3 to 1.0 kpc. The extinction is moderate with A v = 0.24 + 0.06 mag. Low- and high ionization as well as fine-structure lines 
show a complicated kinematic structure probably pointing to a merger in progress. We detect one intervening system at z = 2.18. 
Conclusions. GRB-DLAs have a shallower evolution of metallicity with redshift than QSO absorbers and no evolution in their HI 
column density or ionization fraction. GRB hosts at high redshift seem to continue the trend of the metallicity-luminosity relation 
towards lower metallicities but the sample is still too small to draw a definite conclusion. While the detection of GRBs at z > 4 with 
current satellites is still difficult, they are very important for our understanding of the early epochs of star- and galaxy-formation in 
the Universe. 

Key words. GRBs, individual: GRB 100219A, galaxies: high redshift, galaxies: ISM 



1. Introduction 

Determining the properties of the interstellar medium (ISM) 
in high redshift galaxies is difficult due to their extreme faint- 
ness. Detailed and resolved studies using emission lines in the 
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optical/near-IR range is restricte d to redshifts up to z ~ 3 
with current facilities (see, e.g., |Forster-Schreiber et al.||2009} 
Maiolino et al.| [2008 ) and only with the help of gravitatio nal 



lensing beyond this redshift (see, e.g., Swinbank et al. 2007 for 
resolved studies of a gravitationally lensed galaxy at z = 4.9). 
The use of GRBs as lighthouses, similar to what has been done 
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since many years with intervening absorption systems towards 
QSOs, is a promising way to study the properties of galaxies 
Z > 4. In contrast to QSOs, GRBs with their higher luminosi- 
ties can probe more central and denser regions of those galaxies. 
Long GRBs are connected to massive star- formation (see, e.g., 
|Woosley & Blo om 2006 for a review of the GRB -supernova con- 
nection), he nce GRBs prob e star-forming places across cosmic 
history (e.g., |Fruchter et al ."120 06). 

GRB afterglows constitute a powerful source for absorption 
line studies as they have smooth, featureless intrinsic spectra on 
which the intervening ISM in the host imprints absorption lines 
from different metals. Optical spectroscopy allows the determi- 
nation of metallicities from z ~ 1.7 to ~ 6 when those metal 
absorption lines get redshifted into the near-infra red (NIR). X- 
sh ooter at the ESQ Very La rge Telescope (VLT) (D'Odor ico et 



|aL| |2006| | Vernet et~aT] |20 1 1 \ with its wavelength coverage from 



3,000 to 24,800 A facilitates detailed investigations of galaxies 
up to redshifts of z~ 11-13, provided GRBs exist at such high 
redshifts. The broad wavelength coverage of X-shooter also al- 
lows to fit an extinction curve and to derive the reddening from 
the shape of the continuum, assuming an intrinsic powerlaw for 
the continuum. In most cases, its medium resolution is also good 
enough to allow the fitting of the absorption lines with Voigt 
profiles, allowing a more accurate determination of the column 
density compared to low resolution methods. 

Previous studies have revealed a diverse picture of GRB host 
metallicities but they suffer from a small sample size, with cur- 
rently about 20 metallicities published in the literature (for the 



largest samples see, e.g., |Fynbo et aL]|2006[ |Prochaska et al. 
[20071 |Savaglio et aT]|2009| ). Metallicities of GRB hosts seem 
to be on average higher than those o btained from QSO-DLAs, 
in particular at higher redshifts (e.g. Fyn bo et aT| |2006). GRB 
host metallicities also seem to show a slower evolution with red- 
shift, albeit with a large scatter at any given redshift. Both GRB 
and QSO samples, however, show subsolar metallicities, in some 
cases down to 1/100 solar or le ss (e.g., GRB 050730, |DElia| 
et~aLl|2007l and GRB 090926A, |Rau et al.||2006[ ID'Elia et al.| 



2010|). The difference between the two samples can be explained 



by QSO absorbers probing the outer, less metal-rich regions of 
galaxies due to their lower luminosity and the sight-lines being 
selected according to the cross-s ection, and by GRB hosts being 
larger than QSO-DLA galaxies ( |Fynbo et aLl|2008] ). 

The gas observed in GRB a bsorption-line studi es is mostly in 
a low ionization state (see, e.g., Fynbo et al.|2009| ) which makes 
it a good representation of the cold gas in the host galaxy, mostly 
unaffected by the GRB itself. The same is the case for the ma- 
terial observed in QSO absorbers; however, one has to consider 
that QSO and GRB-DLAs are likely probing different regions 
in the galaxy. Exceptions of highly ionized material in GRB af- 
terglow spectra might point to a special environment of material 
from the vicinity of the GRB itself (e.g., GRB 090426, [Tn5neet| 
|al.|2011| ) and hence are not representative for the metallicity of 
GRB hosts in general. 

Extinction in GRB hosts is usually low, although somewhat 
higher than i n QSO-DLAs. For a sample of 41 Swift GRB after- 
glow spectra, [Zafar et al.|p0TT] ) found that 90% have A y < 0.65 
mag with an average of 0.24 mag and an SMC extinction law. A 
similar result was found from photometric data of Swift GRBs by 
Kann et aL](|2010|). However, high extinction makes the detection 



of the afterglow at optical wave lengths difficult and hence this 
result are somewhat biased (e.g., |Greiner et al.|2011[ Kriihler et 



|al. 12011b] ). The extinction in GRB hosts also seems to decrease 



to the shape of the extinction curve which particularly affects 
optical observations, as these probe the rest-frame UV at those 
redshifts. Furthermore, dust is produced by AGB stars and SNe 
which, in a young star-forming region at high redshifts, might 
not be present or not have produced enough dust yet (see e.g. 
|Gall et al.|20lTT ). 

GRB 1002 19A was discovered in an image trigger by the 
Swi ft satellite ([Gehrels et al 1 12004] ) on 2010 Feb. 19, 15:15:46 
UT (|Rowlinson et al.||2010a|) and had a duration of T90 = 18.8 ± 



5.0 s (Baum gartner et al7| |2010|) . An X -ray as well as an op- 
tical afterglow (Jakobsson et al. 2010| ) was found. A nearby 
galaxy detected in the SDSS was first assumed to be the host 
galaxy (Bloo m & Nugent 1 12010| ), implying a low redshift for 
this GRB. However, GROND (Gamma-Ray burst Optical/Near- 
Infrared Detector) reported a g'-band dropout and suggested a 
redshift of z ~ 4.5 (Krii hler et al.| |2010| ). We took spectra with 
the second-generation instrument X-shooter at the VLT, reveal- 
ing a DLA system and a number of absorption lines at a common 
redshift of z = 4.6667 ( |Groot et al] |2010[ |de Ugarte Postigo 
et al.| 2010a). This redshift was later confirmed by GMOS- 



N/Gemini-N ( |Cenko et al.||2010a| ). 

In the following we present the observations and data reduc- 
tion of the X-shooter spectra (Sect. 2), the analysis of the ab- 
sorption lines in the host as well as the metallicities and abun- 
dances derived from fitting the column densities of the absorp- 
tion lines (Sect. 3), a broad-band fit to the spectral energy distri- 
bution (SED) from X-ray to NIR (Sect. 4) and a description of 
the intervening system at z = 2.18 (Sect. 5). Finally, we place the 
metallicity of GRB 1002 19A in the context of other GRB host 
metallicities and high-redshift galaxies in general (Sect. 6). For 
all calculations we use a cosmology with Q M = 0-27, Q\ = 0.73 
and Hq = 71 km s _1 Mpc _1 . The flux of a GRB is defined to 
behave as F v oc t~ a v~P 



2. Observations 

2.1. X-shooter spectroscopy 

We observed the afterglow of GRB 100219A starting on 2010 
Feb. 20, 02:31 UT (12.5h after the trigger) with X-shooter at the 
VLT under good seeing conditions (-0.7 arcsec). The observa- 
tions consist of 4 different exposures with an exposure time of 
1200 s each. The individual exposures were taken by nodding 
along the slit with an offset of 2" between exposures, using a 
standard ABBA sequenc^] The spectra cover the wavelength 
range from ~ 3000 to 25,000 A at a resolution of R~ 9,000- 
12,000 (visible/NIR range) and an average S/N of ~ 4 per pixel. 
Due to the high redshift and the faintness of the source, we de- 
tect a continuum from the afterglow in the range from -5100 
(Ly break) to 23,000 A. 

The spectra were processed with version 1.3.7 of the X- 
shooter data reduction pipeline ( Goldoni et al. , 2006; Modigliani 
et al. 2010| ). The pipeline transforms the detected counts in 
physical quantities while propagating the errors and bad pixels 
consistently during the process. The reduction consists of the fol- 
lowing steps. First the raw frames are subtracted from each other 
and cosmic ray hits detected and corrected using the method de- 

The frames were then 



veloped by |van Dokkum et al7| (2001 



divided by a master flat field produced using daytime flat field 
exposures with halogen lamps. From the flat-fielded raw frames, 
the orders are extracted and rectified in wavelength space us- 



with redshift (Ka nn et al.J |2010[ |Zafar et al.| |2011| ). Again, we 
might miss more highly extinguished high-redshift bursts, due 



1 The spectrum is taken at a pos. A in the slit, then shifted along the 
slit to pos. B etc. 
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ing a wavelength solution previously obtained from calibration 
frames. The resulting rectified orders are then shifted by the off- 
set used in the observation and merged to obtain the final 2D 
spectrum, in overlapping regions, the orders are weighted by the 
errors. 

From the resulting merged 2D spectrum we extracted a ID 
spectrum together with a corresponding error spectrum and bad 
pixel map using a simple PSF weighting scheme under IDL. The 
error spectrum is derived from the noise of the raw frames propa- 
gated through all the reduction steps using the standard formulae 
of error propagation. However, probably due to the presence of 
the nearby galaxy, the absolute flux level of the ID spectrum is 
too low. Indeed in the UVB arm and in the strong atmospheric 
absorption bands in the NIR arm the median flux goes below 
zero by an amount slightly greater than the 1 cr error which is 
unphysical. We therefore applied the same positive offset to the 
flux in the three arms to correct for this effect and put the median 
flux to zero in those regions. 

To perform flux calibration we extracted a spectrum from an 
observation of the flux standard GD71 (Bohli n et al.| |1995] ) in 
staring mode. The reduction procedure is the same as for the 
science object, with the exception that the sky emission lines are 
subtracted according to the Kelson (2003) method. The standard- 
star spectrum was then divided by the flux table of the same star 
obtained from the CALSPEC HST databas^Jto produce the re- 
sponse function. After division by the flux table, we verified that 
the shape of the response function was compatible with the ones 
found in previous reductions by our group. The response is then 
interpolated in the atmospheric absorption bands in VIS and NIR 
and applied to the spectrum of the source. No telluric correction 
was performed. The final spectrum was transformed from air- 
to vacuum wavelengths. An image of the field including the slit 
position and a part of the 2D spectrum is shown in Fig [I] 

2.2. Imaging observations 

The Gamma-Ray burs t Optical/Near-Infrared Detector GROND 
dGreiner et aL| [2008] ) mounted at the MPG/ESO 2.2 m tele- 
scope at La Silla observatory in Chile observed the field of 
GRB 100219A starting February 20, 2010, 0:30 UTC, 9.24 h 
after the BAT trigger (Kruh ler et al.| |2010| ). Further imaging 
was performed on three later epochs with mean times of 36.6 h, 
61.5 h and 110.7 h after To. In addition, we obtained two 
epochs of NIR observations in the JHK S filters with HAWK- 
I/VLT ( [Kissler-Pattig et al.||2008] ) at mean times of 59.6 h and 
107.8 h after the trigger, respectively. We also obtained a single 
epoch observation of the field using the Auxiliary Port Camera 
(ACAM) on the 4.2 m William Herschel Telescope (WHT) on 
La Palma, starting 8.8 h post-burst. We used the z-band filter, 
with a total exposure of 1800 s. 

All photometric data were reduced in a standard manner us- 
ing pyraf/IRAF, similar to the procedure outlined in Kruhler e"t| 



al. 



P008| ). Given the presence of a bright foreground galaxy 3 
arcsec from the afterglow, Point Spread Function (psf)-fitting 
photometry was used for the photometry in all epochs. The 
stacked images were flux calibrated against GROND observa- 
tions of SDSS fields (Abazai jan et al.||2b09| ) taken immediately 
before or after the field of GRB 10021 9A for the optical gV iV , 
and magnitudes of 2MASS field stars ( [Skrutskie et al.[|2006] ) for 
the JHK S filters. All data were corrected for the G alactic fore- 
ground extinction of E(B - V) = 0.076 according to Schlegel et 
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Fig. 1. Field of GRB 100219A in z-band from WHT/ACAM at 
8.8 h after the burst. The slit position across the afterglow and 
the nearby galaxy is indicated. The lower panel shows the 2D 
spectrum with the spectrum of the afterglow (lower trace) and 
the unrelated galaxy in the slit (upper trace) with some emission 
lines indicated. 



Bohlin (2007), http://www.stsci.edu/hst/observatory/cdbs/cal- 



spec.html 



al. (1998]) assuming a total-to- selective extinction of Ry = 3.08. 
The WHT observations were calibrated via the GROND mag- 
nitudes for stars in the field. For a log of the observations and 
photometry see Table [T] the light curve is shown in Fig. [2] 

Fitting the afterglow light curve with all available opti- 
cal/NIR data simultaneously with an achromatic power-law re- 
sults in a temporal index of a = 1.31^ ^ with a x 2 of 36 for 
30 degrees of freedom (d.o.f.). There is a trend of flattening or 
rebrightening observed at late times (t > 50 h) and we there- 
fore exclude the last two epochs from the fit. Mao et al7| ( |2012 ) 
also found late achromatic bumps in the optical/X-ray lightcurve 
during the first day. 

To determine whether the late afterglow epochs might be 
dominated by a (bright) underlying host galaxy, we reobserved 
the field on Dec. 27 2011 in /'-band with OSIRIS/GTC and a 
total exposure time of 12x120 s. Since the GRB position lies 
near the bright foreground galaxy at z = 0.25 also present in 
the X-sho oter spectrum (see above), we use GALFIT ( |Peng et 
al. 2002) to subtract the contribution of this foreground galaxy. 
The subtraction with a Sersic-profile works very well as shown 
in Fig.|3]and we detect a weak object at the position of the GRB 
with i' = 26.7 ± 0.5 (2 <x, AB magnitudes), using the same ref- 
erence stars as for the afterglow photometry. The host galaxy 
candidate can therefore not have contributed significantly at 2-3 
days post burst and the deviation from the powerlaw decay of 
the afterglow might be due to a late rebrightening. 



3. Absorption lines in the GRB system 

We detect a number of resonant absorption lines from the GRB 
host galaxy, including a DLA, as well as one intervening sys- 
tem. The strongest absorption lines from the host system show 
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Table 1. Late time optical and NIR photometry. Times are midtimes of the observations compared to the burst trigger time. 
Observations in brackets are the stack of the three images above. Magnitudes are given in the AB system and corrected for Galactic 
foreground extinction of E(B - V) =0.076. Upper limits are 3<x. 



date 


exp. time 


instrument 


g' 


r' 


V 


z' 


(h) 


(s) 




(AB mag) 


(AB mag) 


(AB mag) 


(AB mag) 


8.812 


1800 


WHT/ACAM 


— 


— 


— 


21.14+0.017 


9.51639 


460 


GROND 


— 


22.55 + 0.15 


21.26 ± 0.09 


21.06 + 0.11 


9.73361 


460 


GROND 


— 


22.75 + 0.17 


21.29 + 0.07 


20.99 + 0.10 


10.9039 


1460 


GROND 


— 


22.88 + 0.05 


21.51 +0.05 


21.32 + 0.06 


11.4039 


1460 


GROND 


— 


22.88 + 0.06 


21.55 + 0.04 


21.42 + 0.05 


11.9033 


1460 


GROND 


— 


22.99 + 0.06 


21.64 + 0.04 


21.44 + 0.07 


(11.4036 


4380 


GROND 


25.65 + 0.37 


22.96 + 0.04 


21.54 + 0.03 


21.34 + 0.04) 


36.5767 


4380 


GROND 


> 25.65 


24.48 + 0.14 


23.28 + 0.11 


23.39 + 0.17 


61.4692 


4380 


GROND 


> 25.64 


24.99 + 0.20 


23.60 + 0.15 


23.65 + 0.18 


110.692 


4380 


GROND 


> 25.39 


25.40 + 0.41 


24.75 + 0.42 


23.95 + 0.32 


ooj days 


1 a /in 
144U 








ZD. / + U.J 










J 


H 


K s 










(AB mag) 


(AB mag) 


(AB mag) 




10.9106 


1200 


GROND 


20.78 + 0.17 


20.24 + 0.16 






11.4106 


1200 


GROND 


20.99 + 0.20 


20.32 + 0.16 






11.9097 


1200 


GROND 


20.75 + 0.16 


20.65 + 0.25 






(11.4100 


3600 


GROND 


20.84 + 0.12 


20.43 + 0.12 


20.24 + 0.26) 




36.5831 


3600 


GROND 


> 22.22 








59.5956 


2280 


HAWK-I 


23.79 + 0.47 


> 22.84 


22.5 + 0.23 




107.760 


4560 


HAWK-I 


24.08 + 0.34 


23.79 + 0.27 


23.06 + 0.24 
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Fig. 2. Multicolor light curve of the late optical observations 
from GROND, HAWK-I and ACAM. Solid/dashed lines are the 
decay slope of a = -1.31 as derived from a common fit to the 
r', i and z! bands at the first four epochs. The following two 
epochs were excluded as there might be some rebrightening of 
the afterglow. The last point is the detection of the host galaxy. 
Magnitudes are in the AB system and corrected for Galactic ex- 
tinction (see values in Table p} . 



at least 5 velocity components spanning a range of ~ 200 km s" 1 . 
We also detect two fine-structure transitions of Sill and one 
from C II, as have also be en detected in a numbe r of other GRB 
spectra (first detected by |Vreeswijk et al.]|2004| ). These transi- 
tions are likely produced by UV pumping of the excited levels 
by the radiation from the GRB itself. However, Si II* has also 
been observed in star- forming galaxies where they could be pro- 
duced byUV radiation from young stars (see e.g. |Christensen 
et al. ( 2011| )). Finestructure lines that were so far only found in 




Fig. 3. Observations of the host galaxy candidate of GRB 
100291A with OSIRIS/GTC, V band. Left panel: Combined 
/'-band observations. Middle panel: Model of the foreground 
galaxy, the center is indicated by a cross. Right panel: Image 
with the residuals after subtracting the galaxy model. A faint ob- 
ject is visible with /' =26.7+0.5 (circle). 



GRB environments such as Ni II* or Fe II* are not detected in 
our spectrum Christensen et al. ( 2011| ); [Vreeswijk et al. (2007 ). 



3.1. Line fitting and properties 

From the absorption lines dete cted in the X-shooter spectrum 
de Ugarte Postigo et al.| ( |2010a| ) had determined a redshift of 
z = 4.6667 + 0.0005, which we refine in the following. We 
downgrade the spectrum to a resolution of ~ 2 A to eliminate 
the velocity structure and fit Gaussians to all well-detected ab- 
sorption lines. We do not include the high-ionization lines such 
as Si IV and C IV as well as the fine- structure lines since those 
have a different velocity structure (see also Sect. |3.3| ). This gives 
us the redshift of the gas with the strongest absorption com- 
ponent, which coincides with the redshift of the main absorp- 
tion component of S II, the weakest absorption line detected. 
Wavelengths have been corrected for the heliocentric velocity 
shift of 5.26 kms" 1 at the time of our spectrum. The final red- 
shift is z = 4.66723 + 0.00037, which we adopt as v = Okms" 1 . 

The equivalent widths (EWs) are also determined from the 
Gaussians fitted to the lines in the smoothed spectrum. The er- 
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16500 17000 17500 18000 

wavelength (A) 



Fig. 4. X-shooter spectrum of 100219A with the detected absorption lines from the host system (red) and the intervening system 
(green) indicated. Horizontal black bars denote the atmospheric absorption bands, the dotted line shows the error spectrum. The 
UVB arm is not shown since it is below the Lyman limit. The spectrum above - 18, 000 A is not plotted due to low S/N and the 
absence of detected lines. 



rors on the EWs are based on the noise spectrum in the 95% 
width interval of the fitted Gaussian. Errors from the Gaussian 
fit itself are not included. In Table |2] we list the detected lines 
with their corresponding wavelength and EWs, values are given 
in restframe. 

From the full-resolution spectrum we determine the column 
densities by fitting Voigt profiles to the different absorption com- 
ponents using the FITLYMAN context in MIDASjj This pro- 
gram fits Voigt profiles with the atomic parameters for each tran- 
sition by convolving them with the resolution of the instrument. 
Wavelengths of the individual components, column density and 
b-parameter (turbulent and thermal) can ideally all be obtained 
by fitting the absorption profile. The thermal b-parameter is set 
to a low value since temperature broadening is not expected to 
play a role for the observed cold/warm interstellar medium. 

Due to the low S/N, limited resolution and blending of dif- 
ferent components, we obtained the turbulent b-parameter of the 
resonant lines from fitting one wing of the different components 
in the element that, in each specific component, showed the least 
blending and/or saturation. The resulting b-parameters are: 11.2, 
9.5, 10.2 and 6.5 kms" 1 for comp. I, II, III and V respectively 
and 6, 10, 10 and 5 kms -1 for the high ionization lines. For 
comp. and IV we take the average b-parameter of 12.6 km s" 1 
(see A. de Ugarte Postigo et al. in prep.), b-parameters for Mg as 
well as the intervening system could be obtained directly from 



FITLYMAN. Note that the Mg absorption lines, as the only lines 
in the IR-arm of the spectrum, are better resolved than the lines 
in t he VIS-arm due to the slightly higher resolution (see Sect. 
2J}. The b-parameters from fitting blended, unresolved lines 



do not have a real physical meaning ([Jenkins et al. [1986) and 



the low b-parameters obtained (after the spectrum had been con- 
volved with the instrumental resolution) confirms that the ab- 
sorption components are actually unresolved. 

The column density of hydrogen was obtained from fitting 
the red wing of the Lya line (see Fig(5j. Ly/3 is also in the range 
of the spectrum but cannot be fit since it is embedded in the Lya 
forest and the Lya forest blanketing is already playing a role at 
z =4.7. For Mg II, we can only fit the /12803A line of the doublet 
since the /12796A line is heavily contaminated by atmospheric 
lines and even for Mg II /12803A we are only able to give fits for 
the three bluemost components. This is the highest redshift for 
which Mg II has been detected in a GRB host and comparable to 
the highest strong (EW > 1 A) Mg II absorbers detected in QSO 
sightlines (e.g., Jian g et al.|2007| ). 

In Table [3] we present the velocities and column densities 
for the different components. Fig. [6] and [7] show the Voigt profile 
fits to one transition from each ion, selected to either have 
the best S/N or being the least saturated transition or the one 
least affected by skylines and/or blending with other lines. The 
components are not at the same velocities for all lines which 
will be discussed in Sect. 13.31 



http://www.eso.org/sci/software/esomidas/ 
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Table 3. Column densities for the 5 velocity components in the absorption lines as shown in Fig. |6jand|7j Transitions in brackets 
are detected but affected by skylines or blended and have not been used to fit the column densities. 










I 




II 




III 




IV 




V 




Ion, Transition (A) 


Av 


logN 


Av 


logN 


Av 


logN 


Av 


logN 


Av 


logN 


Av 


logN 




(km/s) 


(cm- 2 ) 


(km/s) 


(cm" 2 ) 


(km/s) 


(cm" 2 ) 


(km/s) 


(cm" 2 ) 


(km/s) 


(cm" 2 ) 


(km/s) 


(cm" 2 ) 


SII 1250/53/59 


_ 


_ 


+35 


14.42 ± 0.24 





15.16 + 0.18 


_ 


_ 


_ 


_ 






CII 1334 






+35 


14.99 ± 0.37 





13.78 + 0.26 


-28 


14.24 + 0.28 


-67 


13.60 + 0.21 


-95 


13.87 + 0.36 


CII* 1335.6/35.7 






+35 


13.62 + 0.21 





14.81+0.53 


-28 


14.76 ± 0.27 


-67 


13.62 + 0.18 






CIV 1548, 1550 










+ 15 


13.49 + 0.19 


-18 


13.92 + 0.17 


-45 


13.58 + 0.13 


-100 


14.07 + 0.19 


Sill 1304 (1260, 1526) 


+75 


13.79 + 0.33 


+45 


14.44 ± 0.48 


+ 18 


14.59 ± 0.45 


-13 


14.41+0.35 


-38 


13.91+0.18 






Sill* 1264 






+45 


12.38 + 0.32 


+ 18 


13.22 + 0.26 






-38 


13.25 + 0.19 






Si IV 1393, 1402 










+ 15 


12.53 ± 0.42 


-18 


13.62 + 0.34 


-45 


13.68 + 0.18 


-100 


13.44 ± 0.23 


Fell 1608, (2344) 






+45 


14.11+0.15 


+ 18 


14.25 ± 0.22 


-13 


14.20 ± 0.20 


-38 


13.55 + 0.27 






OI 1302 






+45 


15.93 + 0.43 





14.67 ± 0.43 


-28 


16.76 + 0.56 






-95 


14.18 + 0.25 


A1II 1670 


+75 


12.65 + 0.19 


+45 


13.21+0.39 


+7 


13.32 + 0.48 


-13 


13.00 + 0.65 


-38 


12.96 + 0.26 


-95 


12.88 + 0.36 


MgI2852 






+55 


11.47 + 0.49 


+ 16 


12.24 ± 0.36 


-13 


12.35 + 0.34 


-43 


12.54 + 0.38 






Mgll 2803, (2796) 


(+60 


14.99 ± 0.79) 






(+6 


14.74 ± 0.70) 


-45 


13.38 + 0.14 


-77 


13.28 + 0.63 


-111 


12.94 + 0.17 



Table 2. Absorption line list for the GRB and the intervening 
system. All wavelengths are in vacuum, EWs are given in rest- 
frame. The Mg II lines of the GRB host system are severely af- 
fected by atmospheric lines, hence we cannot give exact values. 



line 


^rest 


^obs 


EW 


logN 




(A) 


(A) 


(A) 


(cm" 2 ) 


GRB 










Ly a 


1215.670 


6894 




21.14 + 0.15 


SII 


1250.584 


7090.375 


0.12 + 0.04 


15.23 + 0.05 


SII 


1253.811 


7107.643 


0.14 + 0.04 


15.23 + 0.05 


SII 


1259.519 


7139.867 


0.10 + 0.04 


15.23 + 0.05 


Sill 


1260.422 


7145.059 


0.52 + 0.04 


15.03+0.06 


Si IP 


1264.738 


7170.786 


0.47 + 0.05 


13.57+0.06 


OI 


1302.169 


7381.448 


0.67 + 0.04 


16.82+0.10 


Sill 


1304.370 


7394.606 


0.55 + 0.05 


15.03+0.06 


CII 


1334.532 


7564.966 


0.56 + 0.05 


15.12+0.06 


CIP 


1335.663 


7571.421 


0.45 + 0.05 


15.11+0.04 


Si IV 


1393.755 


7900.031 


0.40 + 0.05 


14.08+0.08 


Si IV 


1402.770 


7950.811 


0.34 + 0.05 


14.08+0.08 


Sill 


1526.708 


8654.931 


0.50 + 0.07 


15.03+0.06 


CIV 


1548.195 


8775.541 


0.38 + 0.08 


14.43+0.04 


CIV 


1550.770 


8789.557 


0.41 + 0.04 


14.43+0.04 


Fell 


1608.451 


9118.284 


0.33 + 0.04 


14.70+0.03 


A1II 


1670.787 


9471.844 


0.79 + 0.07 


13.84+0.05 


Mgll 


2796.352 


(15852.8) 




(> 13.72) 


Mgll 


2803.530 


15893.554 


-1.4 


> 13.72 


Mgl 


2852.965 


16167.345 


0.45 + 0.08 


12.89 + 0.10 


Interv. 










Fell 


2600.173 


8270.620 


0.10 + 0.04 


13.13 + 0.12 


Mgll 


2796.352 


8895.826 


0.90 + 0.08 


13.69 + 0.10 


Mgll 


2803.530 


8918.049 


0.50 + 0.10 


13.69 + 0.10 


Mgl 


2852.965 


9074.719 


0.10 + 0.05 


12.07 + 0.14 



3.2. Metallicity 

The total column densities of different absorption lines allows 
us to determine the metallicity of the galaxy hosting GRB 
1002 19 A. Leaving out Al and Mg which are affected by atmo- 
spheric lines we get the following metallicities: [C/H] = -2.1 ± 
0.1, [O/H] = -1.0 + 0.1, [Si/H] = -1.5 + 0.1, [S/H] = -1.0 + 0.1 
and [Fe/H] = -1.9 + 0.1. For C and Si we added the contribu- 
tion of the low- and high-ionization resonant and fine- structure 
transitions. 

O and S are the elements known to be least affected by dust 
depletion but OI is usually saturated and the column density 
might not be fully reliable. However , for o ur spectrum they both 
show similar metallicities. Jenkins (2009) investigated the de- 




-4000 -2000 2000 4000 
velocity in (km/s) 



Fig. 5. Lya and Lyfi absorption at the redshift of the GRB host. 
Ly/3 was not fitted due to the Lya forest, the fit shows the Lyf3 
Voigt profile with the parameters adopted from the fit to Lya. 



pletion for a range of different elements in sightlines through 
the MW. With the fits obtained for those elements, they devel- 
oped a method to determine the degree of depletion F* and the 
metallicity by a linear fit of the individual behaviour of the dif- 
ferent elements in case of depletion. The normalization of the 
linear fit then gives a synthetic hydrogen abundance, whose dif- 
ference to the measured H I column density (if available) gives 
the metallicity. Applying this method for our spectra results in a 
mild depletion (F* = 0.15), leaving out the low abundance from 
C, a depletion pattern similar to a warm disk+halo in the MW 
(see also Savage & Sembach 1996 ) and a metallicity of around 
-1.0. In Fig [8] we plot the values from GRB 100219A according 
to the method of Jenkins (2009 ) together with the depletion fac- 
tors of different patterns from our MW as defined in Sava ge & 
|Sembach| ( [l996| ). 
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-200 -100 100 
velocity in (km/s) 



200 



Fig. 6. Absorption components for the resonant transitions in the 
host of GRB 100219A. For Mg II, we show the fit to all com- 
ponents but do not consider the two redmost components for 
the column density due to contamination by the neighboring sky 
line. 



Throughout the paper, we will adopt the metallicity of S 
since it is known to be least affected by dust depletion and there- 
fore assume a metallicity [M/H] = -1.0 ± 0.1 or ~ 1/10 Z© for 
GRB 1002 19 A. This metallicity lies well within the value found 
for other GRB hosts at lower redshifts (see Fig. [9]). In fact, GRB 
host metallicities in the range 3 < z < 4 show a large scat- 
ter between < 1/100 Z Q (GRB 050730,]D'Elia et al.|2007) and 
> 0.3 Z Q (GRB 090205, [D'Avanzo et al|2010| which has only a 
slightly lower redshift than GRB 100219A). QSO-DLA metal- 
licities at this redshift range are considerably lower with a mean 
metallicity of around 1/100 Z© . Measurements from QSO ab- 



x 

M— 

"D 

N 

05 



O 




-200 -100 100 
velocity in (km/s) 



200 



Fig. 7. Comparison between the absorption component of the 
low ionization, fine- structure and high ionization lines. While 
the low ionization lines have the strongest absorption in comp. I, 
it is comp. Ill - IV for the high ionization lines and fine- structure 
transitions. Comp V is rather strong for the high ionization lines 
but is absent in the fine- structure and weaker in the resonant ab- 
sorption lines. 



sorption systems at z = 5 - 6 (Becker et al. 2011 ) indicate an 



even lower metallicity for those galaxies. We will discuss this 
further in Sect. 6.2. 



3.3. Kinematics 

The strongest absorption lines cover a velocity range of around 
150 km s" 1 , but with different components and distributions. 
C II, OI and Al II show the largest velocity spread while Si II 
and Fe II are narrower and have a somewhat different distribu- 
tion. Mg II has a completely different velocity distribution and 
might trace different material than the other resonant lines. 

While for GRBs high and low ionization lines often show 
the same velocity range, the opposite is true for QSO absorbers 
where high and low ionization lines trace different volumes (ion- 
ized halo gas vs. low ionization gas in the central parts of the 
galaxy). The situation is different for GRB 100219A. The high 
ionization lines of Si IV and C IV are blueshifted compared to 
their low ionization counterparts but have a similar distribution 
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Fig. 8. Abundances measured in GRB 100219A and plotted 
according to the dust depletion method described in |Jenkins| 
( 2009 ). A x gives the propensity of an element to convert from 
gas to solid which roughly correlates with element condensation 
temperature and it hence fixed for each element, y shows the 



column density convolved with several parameters that Jenkins 
(2009) obtained by fitting several sightlines in the MW such that 
for a certain depletion along a sightline, the elements lie on a 
straight line with slope F* (the "depletion factor"). The metal- 
licity is then obtained by the difference between log Nh and 
the intersection of the linear fit with the y-axis (the normaliza- 
tion). The solid line shows the linear fit to our data, resulting 
in a depletion factor of F*=0.16 and a metallicity of — 1.0. The 
dashed lines show abundance pattern according to different types 
of sightlines as described in Savage & Sembach|( |1996| ), assum- 
ing the same metallicity of -1.0. C was not included in the fit. 



of components. The mean absorption of the high ionization lines 
is at v = -50kms _1 while it is v = Okms -1 (by definition) for 
the low ionization lines. This also implies a different ionization 
of the material in the individual velocity components with a very 
low ionization for comp I and II, the strongest components of 
the low ionization lines, provided the absorption actually arises 
in the same place for low and high ionization lines. 

The situation is equally puzzling comparing the fine- 
structure lines. The fine- structure lines of Sill* and CII* show 
a similar velocity width than the resonant counterparts but the 
strongest absorption component of the fine- structure line does 
not coincide with the main absorption (or the redmost) com- 
ponent and lies at around -40 kms -1 . Provided that the fine 
structure lines originate from the UV radiation of the GRB, the 
strongest fine- structure absorption should be closest to the GRB 
site and is therefore often found to be in the redmost compo- 
nent (see, e.g., |Thone et al.|2008| ), as expected from an ordered 
velocity field in the host galaxy. If the GRB radiation is indeed 
responsible for the excitation of our fine- structure lines, compo- 
nent III must be the one closest to the GRB as we will show in 
the following. 

Using the ratio between fine- structure and resonant abun- 
dances of the same element, one can determine a rough dis- 
tance of the material from the GRB itself, assuming the lev- 
els are populated by indirect UV pumping from the GRB (see 
Prochaska et al. 2006; Vreeswijk et al.|[2 007). We compute the 
distance between the GRB and the absorbing systems using a 
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Fig. 9. Currently published GRB host metallicities from absorp- 
tion lines (filled dots) including X-shooter metallicities (filled 
squares) compared to QSO-DLA metallicities (diamonds and 
filled triangles), error bars are omitted for clarity. The GRB 
values are taken from Fyn bo et al. ( |2006| ) including lat er pub- 
lished measu rements: GRB 060223 A and GRB 605 10B ( |Chary| 
^FaL| |2007), GRB 060 526 ( |Thone et al.[ |20TTT ), GRB 070802 
dEliasdottir et"all[2009l), GRB 71031, GRB 08 0310 and GRB 
08041 3A JLedoux et al.[ [20091), GRB 080210 ( |De Cia et alT] 
^OTTT), GRB 080607 ( [Prochaska et al.[ [2009] ), GRB 090205 
( [D'Avanzo efall|20l0) , and GRB 090516 (A. de Ugarte Postigo 
in prep.). For GRB 050904 (z = 6.29) we adopt a new, lower 
value, as described in Sect. |6.1[ The values from X-shooter are 



from D' Elia et al-| ( |20T0> (GRB 090926 A), M. Sparre et al. in 



pre p. (GRB 090809, GRB 100425 A) and |de Ugarte Postigo et 
aLl( |2010b| )/J. Fynbo ( priv. comm.) for GRB 09 0313. QSO-DLA 
metallici ties are from |Prochaska et aL] ( [2003 [ ) (diamonds) and 



|Rafelski et aT] ( [20T2| ) (filled triangles) . 



time-dependent, photo -excitation code (see D'El ia et"aLl|2007| 
|Vreeswijk et al. [2007 for details) predicting the ratio between 
fine structure and ground state levels. We use as input the GRB 
flux, spectral and temporal indices estimated in this paper, a 
Doppler parameter of 12.6 km/s (see section 3.1), and we assume 
that the gas is completely in its ground states when the GRB 
explodes. The output of our code is then compared to the Si II 
column densities for the components that also show absorption 
from Si IP, i.e., components I, II and III. The initial Sill ground 
state column densities is obtained by summing the contribution 
from Si II and Si II* for each component. Given the large error 
bars of our data, we only derive an order of magnitude for the 
GRB -absorber distances. We find that component IV is the clos- 
est to the GRB, at a distance of div ~ 300 pc and a similar 
distance for comp. II while component I is at di ~ 1 kpc. 

The complicated velocity structure of the different absorp- 
tion species could be an indication for a chaotic velocity field in 
this galaxy. High-redshift star-forming galaxies are likely noth- 
ing close to the ordered, rotation-supported galaxies we see to- 
day. IFU observations of galaxies at z ~ 2 - 3 show a high veloc- 



ity dispersion of the measured emission lines (see, e.g. Law et 



ST] p009 , and references therein) whose origin is still highly de- 



bated. Among other mechanisms, mergers ([Puech et al.| |2006 
2007) and massive star formation (Lehnert et al., 2009) have 



been proposed. Forster-Schrei ber et al. ( |2009| ) conclude from 
galaxies in the SINS survey that around 30% of the galaxies at 
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z=l-3 are rotationally supported, but turbulent disks probably 
indicate they have completed most of their formation process, 
while another 30% might be galaxy mergers. However, |Law et| 
^T]([2009 ) suggest that clumping of gas during the formation pro- 
cess and subsequent infall to the center of the galaxy can explain 
the observed velocity field and mimic the appearance of a veloc- 
ity fiel d produced by a merger . On the other hand the reverse is 
true as Robertson et al. ( 2006 ) show that a gas-rich merger can 
mimic a velocity field resembling to a disk rotation. Both merg- 
ing and gas infall has to play an even larger role for galaxies at 

The complicated velocity field observed for the host of GRB 
100219A might either be explained with the infall of metal rich 
material onto the galaxy or with the galaxy being involved in a 
large merger with another galaxy. Infalling gas as proposed for 
the velocity structure of other high redshift galaxies should, ac- 
cording to the models, be largely pristine gas, hence this scenario 
seems less likely for our galaxy and the velocity field might in- 
deed point to an ongoing merger. The large velocities observed 
for Mg II could also indicate some sort of galactic outflow, like 
it has been found for GRB hosts at low redshifts (Thon eet al.| 
2007 ). Future studies have to show whether the velocity struc- 



response matrix was used (v.013). We model the observed spec- 
trum with an absorbed power law model. Our fit i ncludes two 



ture of the ISM in high redshift galaxies and GRB hosts is more 
complex than the one of lower redshift galaxies. 



4. SED from X-ray to IR 

The broad wavelength coverage of X-shooter also enables the 
derivation of the absolute extinction curve in combination with 
the X-ray data assuming that the intrinsic aft erglow spectrum is a 
powerlaw or a broken power-law (see e.g. [Watson & Jakobsson| 
|2012| ). The flux calibration of our spectrum can be considered 
reliable between 7000 and 23,000 A, beyond that, the flux of the 
afterglow is very low and we have to disregard this part of the 
spectrum. Below 7000 A, the spectral flux is decreased due to 
the Ly a forest. As a consistency check for the flux calibration of 
the spectrum, we also compare the continuum with the fluxes ob- 
tained from our GROND photometry at the time of the spectrum 
(midtime around 13.2 h after the burst). The spectral continuum 
has been raised by a factor of 1.5 to match the photometric data 
which can be explained by slit losses. As found by our group in 
previous reduction (see e.g. Sparre et al. 2011), the relative flux 
calibration is reliable even at these flux levels while the absolute 
calibration may need to be adjusted with photometry. 

In order to cover a broader wavelength range, we also include 
the X-ray spectrum from XRT on Swift. The XRT began observ- 
ing the field at 15:18:45.0 UT, 178.5 seconds after the BAT trig- 



ger (Rowlinson et al. 2010a). The Swift XRT observed GRB 



100219A from 182 s to 88.5 ks for a total of 8.3 ks. The data are 
mainly in Photon Counting (PC) mode, plus 23 s in Windowed 
Timing mode. The light curve can be modeled with an initial 
power-law decay with an index of a\ = 1.53+Q22* followed by 



a break at 7 85 s to a? = 0.5 7 ± 0.06 and a final steepening at 
~ 4 x 10 4 s Rowlinson et al. ( 2010b| ). A flare is present around 
1300 s and another smaller one at ~ 2 x 10 4 s. 

We extract a spectr um using the Swift XRT spectrum reposi- 
tory ( [Evans et al.] |2009| ) in the time interval 5000-80000 s, dur- 
ing which the X-ray spectral hardness ratio remained constant 
(and therefore we do not expect strong spectral variations), re- 
sulting in 5.7 ks effective exposure time. The resulting X-ray 
spectrum is rebinned to have at least 7 counts per energy bin 
and has been fit using the Churazov weighting scheme within 
the spectral fitting package XSPEC (v. 12.7.0). The latest XRT 



absorber s, one Galactic (fixed to 6.5 x 10 20 cm -2 , |Kalberla et 
al.|20 05 ) and one intrinsic to the GRB host galaxy at z = 4.667 
(we adopt the phabs model). The resulting fit is good with a re- 
duced x^ Qd = 1.1 for 33 degrees of freedom. The photon index is 



r = 



1.57+Q i 3 (90% confidence level for one interesting param- 



eter), the intrinsic column density is Nhx(z) = (3.1 +27) x 10 22 
cm -2 . The mean 0.3-10 keV (unabsorbed) X-ray flux during the 
select interval is 2.5 x 10" 12 erg cm -2 s" 1 . 

The best fit to the X-shooter continuum results in an extinc- 
tion of Ay = 0.24 ± 0.06 mag and an SMC extinction law (see 
Fig [TO]). The moderate extinction from the afterglow continuum 
is in line with the moderate dust depletion derived in Sect. 3.3. 
Dust depletion seems to be common in GRB hosts (see, e.g., 
Savaglio et al.|2006| ) and higher than in QSO-DLAs despite the 
generally low extinction found from SED fits to the afterglow, 
the reason for which is still debated. The slope of the spectral 
energy distribution is /3 opt = 0.60 + 0.12 consis tent with the slop e 
from X-rays of j3 x = 0.57 (0 = T- 1) (see also |Mao et aL| [20121 ). 
To obtain a reasonable value for the electron distribution index /?, 
this slope implies th at the cooling break v c , accor ding to the fire- 
ball model (see e.g., Sari, Piran & Narayan 1998j) has to lie blue- 
wards of the X-shooter spectrum and J3 = (p - l)/2 with p=2.2. 
Fixing the slope in the X-rays to the one from the optical spec- 
trum, we get an intrinsic absorption of N HX (z) = (5.1+2 3) x 10 22 
cm -2 (90% confidence level) and^ 2 ed = 1.1. 

The column density derived from X-rays is about one order 
of magnitude higher than the column density of HI as derived 
from the optical of log N=21.14 (see Table [2]). In fact, the X-ray 
column density is above the average for Swift GRBs but compa- 
rable to the mean for z > 4 GRBs of log Nhx -22.5 following 
the trend of increasing column density with redshift (Cam pana] 



et al. 2009). This is partly an observational effect (at higher red- 



shifts, a larger column density is needed for it to even be de- 
tected) but probably also due to an increased in intervening ab- 
sorbers contributing to the total absorbing column. Gas-to-dust 
and metal-to-dust ratios for GRB 1002 19A are high as it is gen- 
erally the case for GRBs, (Watson & Jakobsson 2012), above 
the values for the local group Zafar et al.||2011 This is usually 
explained by ionization of HI and in particular dust destruction 
by the GRB ( [Watson et aL| |2007l [Campana et al.| |2QQ9| |Fynbo| 
et al. |2009| ). This could also explain the very large difference 
between optical and X-ray column densities (about a factor of 
10), as generally found for GRBs, since X-ray absorption is not 
affected by ionization and hence most metals are found in highly 
ionized states due to the absence of a large column density of 
moderately ionized metals ( [Schady et al.||201 i) . 



5. The intervening system 

Despite the high redshift of the GRB, we only detect one in- 
tervening system at a redshift of z = 2.180 which might be a 
consequence of the low S/N of the spectrum. A second interven- 
ing system at z = 1.8 that had been claimed in de Ugarte Postigo 
|et al.| ( [201 0a| ) could not be confirmed in our refined analysis. 
We are not able to detect any possible absorption lines from the 
nearby galaxy at z = 0.25. The redmost absorption transitions 
from the ISM commonly observed are Ca II XX 3933, 3969 A, 
which, at z - 0.25 still lie within the Lya forest of the GRB ab- 
sorption system, the weaker NaD doublet XX 5890,5895 Ais in a 
part of the spectrum affected by skylines. 



9 



C. C. Thone et al.: GRB 100219A with X-shooter 



C/D 
C 
CD 
"O 

X 
13 



10" 



10" 



10" 7 r 
10 1 



!- 



10 1 



10 16 10 17 
frequency (Hz) 




Fig. 10. Broad-band SED fitting of the X-shooter optical-NIR 
and X-ray spectrum. The black line shows the X-shooter spec- 
trum corrected for an extinction of Ay = 0.24 ± 0.06 mag using 
an SMC extinction law. Orange filled symbols overplotted on the 
optical/IR spectrum are photometric points from GROND (see 
Table [I]). The spectral slope with errors from the X-ray data is 
indicated by the shaded region. The extrapolation of the optical 
slope is consistent with the X-ray spectrum (dotted line). X-rays 
and optical data have been scaled to the corresponding flux at 
the time of the observations of the X-shooter spectrum using in- 
terpolations of photometric data from the X-ray and optical light 
curves. 



The refined redshift for the intervening system according to 
the method described in Sect.[3J]is Zinterv = 2.18116 ± 0.00034. 
For the z = 2.181 system, we securely detect the Mgll doublet 
XX 2796, 2803, Mgl and Fell ,12586, whereas the red part of 
the Fe II /12600 transition is affected by a skyline. Mg II shows 
two components separated by 80 kms" 1 while Mgl and Fell 
only show the stronger of the two components, likely due to the 
lower line strength of those transitions (see Fig. 1 1 and Table[4]). 

The reds hift path covered by our spectrum is of A z ~ 2.6. In 
the past (e.g., |Prochter et al.|2005||Vergani et al.|2009| ) 



an excess 

of both weak and strong intervening absorbers have been found 
along the line of sight of GRBs compared to those of QSOs. 
These works consider an upper redshift limit of z = 2.3 in order 
to compare the GRB statistics with the QSO one, based on the 
SDSS spectra. If we consider this upper limit, the redshift path 
covered by the GRB100219A spectrum is Az = 0.7 and the de- 
tection of one strong Mg II system within this Az is in line with 
an excess of the number density of these systems along GRB 
lines of sight. 

|Sud ilovsky et al. ( 2007) performed also a statistical study on 
the GRB foreground CIV absorbers showing that their number is 
consistent with the one in QSO sightlines. The redshift path cov- 
ered by our spectrum for the CIV absorbers is Az ~ 1.15, while 
the limit on the CIV column density is log N (cm~ ^) -14 and 16 
for b-para meters of 20 and 5 km s" 1 , respectively. Sudi lovsky et| 
al. (2007) found 7 CIV systems in agreement with these limits 
within Az ~ 2.25, whereas we find none. We note that our red- 
shift ra nge for CIV absorbers o f 3.5 < z < 4.7 is different from 
that of [Sudilovsky et al.| ( [2007] ) (1.5 < z < 3) and no statistics 
are available for GRB CIV absorbers within this range. Thanks 
to the large spectral coverage of X-shooter, it will be possible to 
build a statistically significant and homogeneous sample of in- 



tervening absorbers and extend the studies to much higher red- 
shifts. 



Table 4. Column densities for the two velocity components in 
the intervening system at z = 2.181. 





I 




II 




Ion, Transition 


Av 


logN 


Av 


logN 


A 


(km/s) 


(cm" 2 ) 


(km/s) 


(cm" 2 ) 


Mgll 2796, 2803 


+52 


13.10 + 0.08 


-32 


13.56 + 0.10 


Mgl 2852 






-32 


12.07 + 0.14 


Fell 2600 






-32 


13.13 + 0.12 



1.5 
1.0 
0.5 

0.0 
1.5 




x 

(D 
N 

05 



0.5 

0.0 
1.5 
1.0 
0.5 

0.0 
1.5 
1.0 
0.5 
0.0 





Fell(2600) 




Mgl(2852) 



Mgll(2796) 



-300 -200 -100 100 200 300 
velocity in (km/s) 

Fig. 11. Fits to the two velocity components in the intervening 
system at z = 2.181. 



6. GRB 100219A in the context of high redshift 
galaxy abundances 

6.1. GRBs at z> 4 

To date, 260 redshifts of GRBs have been determined (of these, 
198 were detected by Swiftf\ This number includes secure pho- 
tometric redshifts and redshifts determined using emission lines 
from the host galaxy, but not limits from optical/UV detections. 
Only 22 (8.5%) of those 260 bursts are above z = 4 (19 bursts 
or 9.6% for Swift), of which 4 were obtained by photometric 
redshift methods. 8 (3.1% or 4.0% for Swift only) are at z > 5 
of which 2 are photometric redshifts. The current records are 
z = 8.2 for GRB 090423A (spectroscopic, jSalvaterra et al.|2009| 



4 http://www.mpe.mpg.de/ jcg/grbgen.html and the GCN archive as 
of 11 June 2012 
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Tanvir et al.|2009) and z = 9.4 for GRB 090429A (photomet- 



ric, |Cucchiara et al.|201 1| ), which makes them one of the farthest 
objects detected in the Universe. 

In the past, it had been estimated that a much larger fraction 
of GRBs might be observed at redshifts beyond 4-5 ( |Bromm| H 
& Loeb| |2002| |2006|) and, taking into account the sensitivity of 



Swift, at least 10% should be at z> 5. It had been suggested that 
part of the high-redshift GRBs were missed due to the lack of 
NIR observations and high sensiti vity NIR spec trographs to de- 
termine redshifts. However, Greiner et al. (2011) detected after- 
glows with GROND for more than 90% of the GRBs if obser- 
vations were possible early enough, although they did not detect 
e.g. GRB 090429B at z=9.4 ( [Cucchiara efS!H25TT) . Their red- 
shift distribution is similar to the overall distribution, leaving lit- 
tle room for a significant high-redshift population being missed 
by observations. X-shooter has to date determined 26 redshifts 
(2 or 8% at z>4) and the redshift distribution for GRBs observed 
with X-shooter is similar to the distribution of all GRBs (see 
Fig. [Hf Observational constraints would also not explain the 
scarcity of GRBs between redshift 4-6, which are still easily 
detectable with optical spectrographs. 

The reason for this low rate might be related to the capability 
of Swift (and other past and current y-ray observatories) to de- 
tect high-z GRBs. Time dilation stretches the emitted flux over a 
longer time interval, making rate triggers more difficult (in fact 
GRB 100219A was discovered in an image trigger). In addition, 
the peak emission of long, already intrinsically soft GRBs gets 
shifted to even lower ba nds, so that part of the flux falls below 
the BAT bands. Infact, |Gorosabel et al.| ( |2004| concluded that 
the INTEGRAL satellite with it's different energy bands should 
detect more GRBs at high z than Swfit. Several authors have also 
suggested that without a luminosity evolution of GRBs it is un- 
likely that we would detect GRBs at that redshift at all (e.g., 
|Salvaterra et al.|2007] ), but determining the luminosity function 
of high redshift GRBs suffers from small number statistics. 

Another unknown factor is the evolution of the star- 
formation rate (SFR) at high redshifts. The peak of the SF is 
expected to be around z ~ 2, and it is not yet known how the 
SF evolves at redshifts beyond z ~ 4. Several attempts have 
been made to determine the cosmic S FR using GRBs (see, e.g., 
Yiiksel et al.|2008|[Kistler et al.|2009] ) indicating a shallower de- 



cay towards higher redshift than derived by galaxy SF measure- 
ments. The redshift distribution of Swift GRBs at high redshifts 
can be well explained by a cons tant cosmic star-formati on den- 



sity from redshift 3.5 onwards ( Jakobsson et al.[ 2006) ^ This 



could be explained if GRBs preferentially select low-metallicity 
regions. Which of the effects, an evolution of the luminosity 
function, a metallicity selection effect or a larger SFR at higher 
redshifts, is the one responsible for the observed evolution can- 
not be determined at present. 

In Table [5] we list the properties of all GRBs with z > 4 
with spectroscopic observations since the discovery of GRB af- 
terglows. Of the 16 GRBs at z > 4 with spectroscopic redshifts, 
5 have measured metallicities, 3 more have derived limits. For 
GRBs with 4 < z < 5, the majority has values or limits for the 
metallicity (7 out of 11), for the remaining 4 GRBs at least two 
seem to have good enough spectra to determine a metallicity but 
are not yet publ ished (GRB 100302A , |Chornock et al.|2010| and 
GRB 100513A, |Cenko et al.||2010b| ). The data of GRB 090516 



5 Since the first submission of this publication, X-shooter had de- 
tected a GRB with higher redshift th an GRB 1002 19 A, GRB 11 1008 A 
at z=4.9898 ( [Wiersema et aLl [20TT] ) 
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Fig. 12. Redshift distribution of all GRBs with secure redshifts 
(spectroscopic and photometric, orange) and the same for GRBs 
observed with X-shooter (blue histogram). X-shooter spectra 
where the redshift could not be determined are shown in red. 



which are available to our collaboration will be published in a 
forthcoming paper (A. de Ugarte Postigo et al., in prep.) and 
we list the metallicity determined from Sill A 1808. This spec- 
trum shows unusually strong lines, so that even the weak Si II 
1808 line is mildly saturated, so the real metallicity is proba- 
bly ~ 0.3-0.4 dex higher, assuming typical values for the b- 
parameter around 20-30 km s" 1 . 

At z > 5, determining metallicities gets increasingly difficult, 
only one metallicity has been reported for a GRB beyond z = 5, 
GRB 050904 ( [Kawai et al.l |2006). However, they list metallic- 
ities derived from C, Si, O and S that differ by more than an 
order of magnitude and they adopted the value from S which 
is the highest from all 4 elements. They argue that this might 
be due to dust depletion, however, e.g., Si should track S rather 
close. |Kann et al.| ( |2007| ) and |Zafar et al.| (|2010), however, found 
no indication for dust depletion from the SED of the afterglow 
itself (but see |Stratta et al.|2007||20TTT ). 

We reinvestigate the spectra of GRB 050 904 from|Kawai et| 
|aL] ( |2006| ) which have also been studied in |Zafar et al.| ( |2010| ) 
and find different EWs for S II A 1253 and Sill A 1260 (2.7 and 
13.0 A observer frame respectively), while S II A 1259 is blended 
with Si II A 1260 an d cannot be measured. All lines except S II 
listed in |Kawai et~aL| ( |2006| ) are completely saturated, also in- 
cluding another Sill line at Al 304 not listed in the paper. The 



saturation makes a curve-of-growth fit as mentioned in Kawai et 



al. (2006 ) impossible, however, the column densities of Si II, O I 



and C II are definitely higher than listed in the paper. Taking the 
EW from S II A 1253 and assuming a (resonable) Z?-parameter of 
~ 35 km s" 1 , we get a column density of aro u nd log N ~ 15.3, 0.3 
dex lower than listed in |Kawai et"aL] ( |2006| ). |Totani et al.] ( |2006| ) 
reinvestigated the DLA column density of GRB 050904 and find 
a slightly higher value of log N#/=21.6. 

Together, this gives a revised metallicity of [M/H] — 1.6 and 
we give a conservative 3cr error of 0.3 dex. This is also in line 
with the findings of |Campana et al.|(2007| ) deriving a metallicity 
of [M/H]>-1.7 from X-ray spectra. Dust depletion is likely not 
playing a major role since the the metallicities from different 
elements are more similar than previously thought, although the 



For a recent update see http://www.raunvis.hi.is/~pja/GRBsample.htmldefinite values cannot be determined. 
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6.2. Comparison between high redshift GRBs and QSO 
absorbers 

Observations of QSO absorbers at z > 4 become difficult since 
current surveys only detect the most luminous ones. The inter- 
galactic medium at those redshifts is almost opaque for Lymana 
photons. This makes it difficult to detect and measure the col- 
umn density of DLAs in the sightline, which in turn hampers the 
direct determination of the met allicity. The largest samples of 
QSO-DLAs and th eir properties ([Ledoux et al . , 2006 ; Prochaska 



eTaL] [20031 [2005 



2006 ) comprise 



Dessauges-Zavadsky et al. 
a number of systems up to z ~ 4.5. Becker et al. ( 2011| ) deter- 
mined properties of 5 < z < 6 QSO absorbers but could only 
indirectly derive the metallicity of those systems. 

QSO-DLAs show a decrease in metallicity with redshift 



(Prochaska et al. 2007 1 Dessauge s-Zavadsky et al.|[2006| ) with 
a large scatter at each redshift. There seems to be a "metallicity 

floor" at ~ [Z/H] 2.6 explained by a rapid early enrichment 

due to star-formation. In Fig. [9] we plot all GRB host metallici- 
ties reported in the literature together with metallicities from the 
large QSO-DLA sample of |Prochaska et al.|p003] ) and fit a lin- 
ear evolution to th e data for eac h sample. As it has been pointed 
out in the past ( [Savaglio efaTl |2006| [Pynbo et al.| |2006), GRB 
hosts have larger metallicities and less evolution with redshift 
compared to QSO-DLAs. At z > 4 this becomes clearly evident 
as the mean GRB-DLA metallicity is about one order of mag- 
nitude larger than the mean of QSO-DLAs. Albeit a large scat- 
ter exists at all redshifts, the mean GRB -host metallicity seems 
to show a much shallower evolution compared to that of QSO- 
DLAs. 

At z > 5, the ionization fraction of QSO-absobers very low 



and Si IV and C IV are hardly ever detected (Becker et al. 2011 



However, this is consistent with the ionization fraction at lower 
redshifts due to a general decrease in column densities which 
makes the column densities of Si IV and CIV fall below the 
detection limit. GRBs with their generally higher metallicity 
should therefore show more Si IV and C IV, which is the case. 
GRBs do not exhibit any trend with redshift, neither in the over- 
all abundance, nor is there any evident change in the ionization 
fraction and at high redshifts as shown in Fig. 



13 for carbon. 



Special cases such as the highly ionized material around GRB 
090426 (not shown in the plot) do not reflect this global picture 



since they likely probe very special environments (Tho ne et al. 

|20TTT ). 

The distribution of GRB-DLA HI column densities have a 
higher mean value tha n those of QSO-DL As and extend to much 
higher densities (e.g. |Fynbo et al.|[2009| ). It has been suggested 
that the cosmic neutral hydrogen density reaches a maxim um 
around z ~ 3.5 and decreases for lower redshifts (Prochaska et 



al. 2005 ), since the gas is consumed or expelled from galaxies, 
while the evolution at z > 4 is still inconclusive. The frequency 
distribution of DLAs and sub-DLAs shows little evolution but 
the slope might change slightly implying more sub-DLAs and 
less high column density systems at high redshifts (Prochaska et 
aTJ [20U5| |Guimaraes eTaL] [2009] ). In Fig. [14] we plot log N H i 



for all GRBs for which this value is available from the literature 
and determine the mean and standard deviation at three redshift 
bins excluding sub-DLAs (log N#/ < 20.3). The mean column 
densities are log N HI = 21.76 ± 0.38 for 2 < z < 3, 21.80 ± 0.62 
for 3 < z < 4 and 21.42 ± 0.56 for z > 4, hence there is no 
evolution with redshift. Only a few GRBs with sub-DLA column 
densities have been found, all at z < 4, but this is consistent 
with the number expected from the overall HI distribution. It 
seems GRB hosts have a high content of neutral gas irrespective 



• z > 4 
■ 3<z<4 
a 2<z<3 



of redshift which gives them a large reservoir for the formation 
of stars. 
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Fig. 13. Ionization of C in GRB -DLAs at 3 different redshift bins 
(filled symbols) and the comparison to measurements in z = 5-6 
QSO absorbers (empty circles) from Becker et al. ( 2011|). GRB- 
E Ws are t aken from the low resolution Swift sample of |Fynbo et| 
[aL] ( |2Q09| ) and its extension until Sept. 2009 published in A. de 
Ugarte Postigo et al. in prep. The CIV EWs for the GRB sample 
include both lines of the AA 1548,1550 doublet, upper limits of 
the QSO sample are 2cr limits on CIV A1548. 
Values for QSOs for the othe r two redshift bins ar e not shown for clar- 
ity and we refer to Fig. 13 in |Becker et al.| ( p011| ). GRB -DLAs do not 
show any evident trend with redshift, neither in line strength nor in the 
ionization fraction and have higher line strengths than QSO absorbers 
at z >4. 



It is still debated whether QSO-DLAs (and GRB -DLAs) 
show an enhancement in a elements (O, Ne, Mg, Si, Ca), pro- 
duced in massive stars, versus iron-peak elements (Fe, Ti, Zn), 
which primarily come from Type la SNe. An overabundance 
of a elements is expected for young star-forming regions and 
a corresponding correlation between metallicity and [a/Fe] has 
been found for MW stars. Determining the [a/Fe] ratio is com- 
plicated by differential dust-depletion since Fe is more sensitive 
to dust depletion than e.g. Si. All studies of QSO-DLAs have 
found a statistically significant a enh ancement (e.g., Wolfe et al. 



2005 ; Dessauges-Zavadsky et al. 2006 ). It has been claimed that, 
correcting for dust depletion, the enhancement disa ppears (e.g., 
Pett ini et al.|2 002 but see also Savaglio et al.|2000| ). However, a 
mean positive ratio of [a/Fe] was found even for very metal-poor 
(hence dust- free) systems ( Wolf e et al. 2005|) as well as for dust- 
free syst ems as determined by [Fe/Zn] ([Dessauges-Zavads ky et| 
|al.|20 06, Zn is supposed to be little affected by dust depletion). 
Extinction determined from the spectral continuum is gener- 
ally low for GRB-DLA and QSO-DLAs with somewhat higher 
values for GRB -DLAs. However, dust depletion and extinction 
are not always tightly correlated, which could be caused by 
a "patchy" distribution of the dust along the line-of- sight. For 
GRB 100219A, we get a ratio [or/Fe] of +0.39 (Si) and +0.76 
(S). Si is only mildly depleted while S can be considered as little 
affected by depletion so the enhancement might be real. An en- 
hancement has also been found for other GRBs ( Savaglio et al. 



2006 ) but the sample is still too small to make a definite conclu- 
sion or to investigate a possible redshift dependence. 
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Fig. 14. Colu mn density of Swift G RB -PL As vs. redshift from 
the sample of |Fynbo et al. (2009 ) and later published values 
from the literature. The vertical dotted line indicates the detec- 
tion limit of Ly a from the ground. Grey points are column den- 
sities, overplotted are the average column density and standard 
deviation in 3 redshift bins (2<z<3,3<z<4 and z > 4). For 
the average column density we exclude sub-DLA systems (log 
Nhi < 20.3, horizontal dotted line) since they might probe dif- 
ferent environments; the values are plotted with empty circles. 
The GRB-DLA HI column density shows a large scatter at all 
redshifts but no trend towards lower column densities at higher 
redshifts as observed for QSO-DLAs. 



6.3. The mass-metallicity relation at high redshifts 

In the local Universe, galaxies follow a correlation between stel- 
lar mass content and metallicity and, with a less tight correlation, 
between (#-band) luminosity and metallicity (L-Z and M-Z re- 
lations). A corresponding mass-metallic ity relation has been de - 
rived from the large SDSS sample by |Tremonti et al.| ([2004), 
with metallicities derived from strong emission lines of the lu- 
minous gas in the galaxy. Deriving a reliable stellar mass content 
of galaxies usually requires broad-band SED fitting or at least a 
detection in the rest frame K band, which makes it difficult to 
derive the masses at higher redshifts. Emission-line metallicity 
measurements at redshifts beyond z ~ 1 require the use of high- 
sensitivity NIR spectrographs. In the past years, the local M-Z 
(or L-Z) relation has been extended to redshifts of z ~ 3.5 (see 
Mannucci et al. 2009 and references therein). 

It has become clear that the M-Z relation at increasing red- 
shift changes towards lower metallicities for a given stellar mass 



([Kewley & Ellison, 2008; Sava glio et all|2005l[ErFeTaLl |2006; 
Mannucci et al.||2009| ). This effect is largely due to the enrich- 
ment with metals over time but the metal enrichment of galax- 
ies is also a complex function of star-formation efficiency and 
feedback with the intergalactic medium. GRB hosts fall below 
the mass-metallicity relation for field (SDSS) galaxies ( |Han et 
0] |2010l |Kocevski & West| [20TT| ) which had been taken as 
evidence for GRB hosts being biased towards lower metallici- 
ties, possibly because GRBs need low-metallicity environments 
( [Woosley & Heger|[2006| ). However, |Mannucci et al.| ( [2010| ) es- 
tablished a connection between metallicity and star-formation 
rate showing that star-forming galaxies have on average lower 
metallici ties, which can explain p art of the discrepancy to field 
galaxies ([Kocevski & West| [201 1 ). Similar to the M-Z relation, 



the L-Z relation also shows an evolution with redshift, but the 
overall correlation is less tight. At low redshifts, the L-Z cor- 
relation is better constrained, probably due to a tight correla- 
tion between mass and luminosity. At higher redshift, the scatter 
in the mass-to-light ratio (M/L) between #-band magnitude and 
mass is rather large (see Er b et al.|2006| ) and also steeper than in 
the local Universe, which makes the L-Z relation at z > 2 rather 
unconstrained. 

Metallicities of GRB hosts using emission lines have only 
been determined for a few hosts up to z ~ 2, which will be sig- 
nificantly improved with X-shooter in the near future. At high 
redshifts, we have to rely on absorption-line metallicities, but 
we have to be aware that the metallicity derived from the lumi- 
nous, hot gas of star-forming regions might not be the same as 
derived from the cold gas of the galaxy ISM. We collected all re- 
ported host detections and upper limits for z > 4 GRBs from the 
literature (see Table[5]). Due to the lack of broad-band modeling, 
we only investigate the L-Z relation since the relation between 
5-band magnitude and stellar mass seems to change with red- 
shift and is still uncertain at higher redshifts, therefore deriving a 
stellar mass would introduce additional errors. For GRB 050904, 
GRB 060223 A, GRB 0605 10B an d GRB 090423, rest-fra me B- 
band magnitudes were given in C hary et aL] ( |2007| |2009j ) from 
Spitzer observation s. Fo r GRB 1002 19 A, we take the /'-band de- 
tection from Sect. |2.2| Magnitudes given in other bands were 
transformed to rest- frame B-band magnitudes using a rough k- 
correction of the form k = -2.5(1 + J3) x /ogio(l + z(^bI knit)) 
with Ainit. and A B as the central wavelength of the original and 
the B-band filter respectively. We assume a powerlaw slope for 
the host SED of the form F v oc with = 0.5. 

Fig. [15] shows the result for GRB hosts with z > 4 com- 
pared to two other sample s for field galaxies of z ~ 2.3 (Erb 
|et~aLl [20061 ) and z ~ 3.1 ( [Mannucci et all [20091 ). GRB 050904 
at z = 6.29 with our corrected metallicity and GRB 060223A 
clearly lie below the galaxies at z ~ 3.1 while the other four 
hosts, including GRB 1002 19 A, are consistent with the relation 
at lower redshifts. The number of high-redshift GRBs is still 
too low to make any conclusive statement on the L-Z relation 
of GRB hosts. 
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Fig. 15. Luminosity-metallicity relation of z > 4 GRBs com- 
pared to field galaxy samples at slightly lower reds hifts. Galaxy 
dat a are fr om |Erb et al.| ( |2006| ) for z ~ 2.3 and from Mannu ccTet| 
|aT1( |2009l )forz~ 3.1. 
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Salvaterr a et al.| ( |2011| ) studied the redshift evolution of the 
mass-metallicity relation from redshift 5 to 10 using galaxy evo- 
lution simulations with the GADGET code ( |Springel| 2005 ). 
They found considerable enrichment even of the smallest galax- 
ies of up to 1/10 already at redshifts 7-8 and probably even 
higher. The metallicity of ~ 0.1 Z for GRB 100219A is fully 
consistent with chemical evolution models. Comparing the L-Z 
relation at z = 5 from the simulations of Salvaterra et al.| ( |20TT] ), 
our GRB hosts at z = 4 - 5 would still fall somewhat below the 
relation, confirming the trend seen at lower redshifts. Further in- 
vestigations of the relation with new observations, in particular 
at z> 5 would be highly desirable to study the early chemical 
enrichment of galaxies and their star-formation history and ISM 
feedback. 

A direct determination of the L-Z/M-Z relation of QSO ab- 
sorbers is difficult to obtain, however, several authors have con- 
firmed a correlation of the velocity spread of, e.g., Sill X 1526 
and the metallicity ( |Ledoux et aLl|2006[|PTochaska et alT]|2Q08j ). 
This correlation seems to decrease with redshift, mostly due to 
the change in metallicity. The velocity spread is probably re- 
lat ed to t he dark-matter halo mass of the absorbers ( [Ledoux et| 
aL||2006), hence the velocity-metallicity relation can be seen as 
a mass-metallicity relation and the observed trend fits to the M- 
Z/L-Z relations found in emission. A tight correlation has also 
been found for the metallicity vs. the EW of Sill X 1526 of 
QSO-DLAs such that it can be used as a prox y for metallicity 
if no other measurement is available (Prochas ka et al.[ |2008| ) 
and this correlation might also be present for GRB-DLAs, al- 
beit with a different slope and offset. The reason for this corre- 
lation might be the underlying velocity-metallicity relation and 
the suggestion that the EW of saturated absorption lines traces 
the kinematics of the gas well. The EW of Sill X 1526 of GRB - 
DLAs does not evolve much with redshift (Fynbo et al. 2009 ), 
while it decreases for QSO-DLAs. On the other hand, the ve- 
locity spread of, e.g., Sill X 1526 in QSO-DLAs does not show 



lar ge cha nges at high redshifts ( [Prochaska et al.| |2005 ; Bec ker et 
|aL||2011| ). The identification of the galaxy counterparts of QSO- 
DLAs is difficult due to their faintness and the proximity of the 
QSO, therefore, a direct mass- or luminosity-metallicity relation 
for the galaxy counterparts of QSO-DLAs has not been deter- 
mined yet. 

7. Conclusions 

In this paper we presented X-shooter spectra of GRB 100219A at 
Z - 4.667 which is the highest redshift at which a high-resolution 
spectrum has been available until now, allowing a detailed anal- 
ysis of the abundances in a high-redshift (star-forming) galaxy. 
The metallicity as determined from S is moderately high with 
[M/H] = -1.0 or 0.1 Z©, similar to the value found for other 
GRBs between redshift 3-4, but 10 times higher than the aver- 
age metallicity of galaxies found in the sightlines of QSOs at this 
redshift. There is little evidence for extinction from the afterglow 
SED and the relative abundances of the detected metal absorp- 
tion species indicate only a mild dust depletion. This suggests 
that the a-element enhancement of [a/Fe] =0.3-0.7 we find from 
Si and S is likely real. The kinematics of low-, high-ionization 
and fine-structure lines is rather complicated, probably showing 
an early galaxy in the process of formation or merging with an- 
other galaxy. We also detect one intervening system at z = 2.181, 
consistent with the detection rate of intervening absorbers in 
GRB sightlines. 

Studying galaxies and their abundances in the high-redshift 
Universe is of great importance to our understanding of the cos- 



mic chemical evolution, in particular at redshifts where galax- 
ies were still in the process of formation. Over the last years it 
has become evident that galaxies probed by QSOs and GRBs 
show somewhat different properties. While GRB hosts have a 
very shallow metallicity evolution, QSO-DLAs have on average 
a factor of 10 times less in metallicity at z ~ 4 and show a steep 
increase up to z ~ 2, where the metallicities of both samples 
become similar. GRBs also show no evolution in ionization rate 
(though this is likely also the case for QSO-DLAs) and in the 
HI average column density and distribution. Probably, GRBs se- 
lect similar types of galaxies across the history of the Universe, 
highly star-forming and metal-poor, while QSO-DLAs probe the 
average galaxy population. 

It has been suggested that both QSO- and GRB-DLAs might 
come from the same population of galaxies and the metallic- 
ity distribution can be explained by a sightline effect (GRBs 
probe the dense regions of the galaxy while QSO-DLAs probe 
the more metal-poor outskirts) combined with a slightly higher 
average mass for GRB-DLA galaxies. This seems likely since 
GRB-DLAs also have much stronger lines, hence probing denser 
regions, and a higher velocity width than QSO-DLAs. It is still 
unclear into what type of galaxies in the present Universe high- 
redshift QSO-DLAs develop into, though they are likely not 
dwarf galaxies ( [Prochaska et al.| |2005| ). GRB hosts at low red- 
shifts are a diverse mix of star-forming galaxies, from blue com- 
pact dwarfs to spiral galaxies. 

The picture is somewhat inverted for the mass-metallicity re- 
lation of GRB hosts compared to field galaxies. At low redshifts, 
GRBs select galaxies more metal-poor for their mass than aver- 
age, even after accounting for a possible bias from star-forming 
galaxies (that have shown to be on average more metal-poor). 
The comparison at higher redshifts is still unclear but the dif- 
ference might be less pronounced than in the local Universe. It 
will be an important task to develop a common picture for the 
galaxy populations probed by different methods in order to get 
the global picture of galactic chemical evolution. 

The evolution beyond redshift ~ 5 is still largely unknown 
since the current detection rate of such high-redshift GRBs is 
still low. At redshifts approaching the reionization epoch and 
the formation of the first stars and galaxies, the determination 
of metallicities and relative abundances are crucial to test the 
models for high-redshift star-formation and the formation and 
evolution of galaxies. 
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